Linear energy-preserving integrators for Poisson systems
نویسندگان
چکیده
For Hamiltonian systems with non-canonical structure matrix a new class of numerical integrators is proposed. The methods exactly preserve energy, are invariant with respect to linear transformations, and have arbitrarily high order. Those of optimal order also preserve quadratic Casimir functions. The discussion of the order is based on an interpretation as partitioned Runge–Kutta method with infinitely many stages.
منابع مشابه
Geometric Exponential Integrators
In this paper, we consider exponential integrators for semilinear Poisson systems. Two types of exponential integrators are constructed, one preserves the Poisson structure, and the other preserves energy. Numerical experiments for semilinear Possion systems obtained by semi-discretizing Hamiltonian PDEs are presented. These geometric exponential integrators exhibit better long time stability p...
متن کاملEnergy-preserving integrators for stochastic Poisson systems
A new class of energy-preserving numerical schemes for stochastic Hamiltonian systems with non-canonical structure matrix (in the Stratonovich sense) is proposed. These numerical integrators are of mean-square order one and also preserve quadratic Casimir functions. In the deterministic setting, our schemes reduce to methods proposed in [9] and [6].
متن کاملA derivation of energy-preserving exponentially-fitted integrators for Poisson systems
Exponentially-fitted (EF) methods are special methods for ordinary differential equations that better catch periodic/oscillatory solutions. Such solutions often appear in Hamiltonian systems, and in view of this, symplectic or energy-preserving variants of EF methods have been intensively studied recently. In these studies, the symplectic variants have been further applied to Poisson systems, w...
متن کاملReversible measure-preserving integrators for non-Hamiltonian systems.
We present a systematic method for deriving reversible measure-preserving integrators for non-Hamiltonian systems such as the Nosé-Hoover thermostat and generalized Gaussian moment thermostat (GGMT). Our approach exploits the (non-Poisson) bracket structure underlying the thermostat equations of motion. Numerical implementation for the GGMT system shows that our algorithm accurately conserves t...
متن کاملEnergy-preserving methods for Poisson systems
We present and analyze energy-conserving methods for the numerical inte-gration of IVPs of Poisson type systems that are able to preserve some Casimirs.Their derivation and analysis is done within the framework following the ideasof Hamiltonian BVMs (HBVMs) (see [1] and references therein). The pro-posed methods turn out to be equivalent to those recently derived in [2], giv...
متن کامل